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Abstract 

Introduction: Onychophora is a relatively small phylum within Ecdysozoa, and is considered to be the sister group 
to Arthropoda. Compared to the arthropods, that have radiated into countless divergent forms, the onychophoran 
body plan is overall comparably simple and does not display much in-phylum variation. An important component 
of arthropod morphological diversity consists of variation of tagmosis, i.e. the grouping of segments into functional 
units (tagmata), and this in turn is correlated with differences in expression patterns of the Hox genes. How these 
genes are expressed in the simpler onychophorans, the subject of this paper, would therefore be of interest in 
understanding their subsequent evolution in the arthropods, especially if an argument can be made for the 
onychophoran system broadly reflecting the ancestral state in the arthropods. 

Results: The sequences and embryonic expression patterns of the complete set often Hox genes of an 
onychophoran {Euperipatoides kanangrensis) are described for the first time. We find that they are all expressed in 
characteristic patterns that suggest a function as classical Hox genes. The onychophoran Hox genes obey spatial 
colinearity, and with the exception of Ultrabithorax {Ubx), they all have different and distinct anterior expression 
borders. Notably, Ubx transcripts form a posterior to anterior gradient in the onychophoran trunk. Expression of all 
onychophoran Hox genes extends continuously from their anterior border to the rear end of the embryo. 

Conclusions: The spatial expression pattern of the onychophoran Hox genes may contribute to a combinatorial 
Hox code that is involved in giving each segment its identity. This patterning of segments in the uniform trunk, 
however, apparently predates the evolution of distinct segmental differences in external morphology seen in 
arthropods. The gradient-like expression of Ubx may give posterior segments their specific identity, even though 
they otherwise express the same set of Hox genes. We suggest that the confined domains of Hox gene expression 
seen in arthropods evolved from an ancestral onychophoran-like Hox gene pattern. Reconstruction of the ancestral 
arthropod Hox pattern and comparison with the patterns in the different arthropod classes reveals phylogenetic 
support for Mandibulata and Tetraconata, but not Myriochelata and Atelocerata. 
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Introduction 

Arthropod segmentation - its origin and maintenance - is 
among the key topics of evolutionary developmental (Evo- 
Devo) research. Arthropods are of particular interest as 
they comprise the most speciose and disparate animal 
phylum. They form a clade together with the onychopho- 
rans and tardigrades, referred to as Panarthropoda [1]. 
The arthropods themselves include four or five (depend- 
ing on the status of pygnogonids) classes [2,3]: the insects, 
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the crustaceans, the myriapods and the chelicerates (with 
or without the pyenogonids), although there is significant 
molecular support for insects being an ingroup to a para- 
phyletic "Crustacea", the group in total being referred to 
as the Pancrustacea or Tetraconata [3,4], While the body 
plans of many arthropod groups have diversified signifi- 
cantly since the Cambrian Period, the lobopodian body 
plan, represented today by the onychophorans, appears to 
have changed relatively little [5-8]. Many of the key char- 
acteristics of the arthropods - such as jointed limbs and 
full adult body segmentation with pronounced segmental 
indentations - are not present in onychophorans, and, 
based on reconstructions of the arthropod stem-group, 
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these absences are best interpreted as being primitive [9]. 
However, traces of such features, such as the genetic tool- 
kit required for podomere patterning are present in the 
limbs of onychophorans [10], and arthropod-like append- 
ages are present in stem-group arthropods such as the 
anomalocaridids [11,12]. 

In Drosophila and to some extent in other arthropods, 
segmentation is under control of a hierarchic segmentation 
gene cascade reviewed in [13-15]. This cascade controls the 
expression of the Hox genes, which in turn specify segmen- 
tal identity. It is believed that the Hox genes are involved in 
providing positional information in a combinatorial mode to 
give each segment its identity along the anterior-posterior 
body axis [16-19]. Thus, disturbance of the Hox patterning, 
such as loss-of-function or ectopic expression of Hox genes 
often results in homeotic transformations, the change of 
one segments identity into that of another [20-23]. Beyond 
that, Hox genes are believed to be involved in tagmosis, i.e. 
the grouping of segments into functional units (tagmata) 
e.g. [24,25]. The Drosophila Hox clusters contain eight Hox 
genes, but the ancestral arthropod Hox cluster most likely 
contained ten Hox genes [26,27]. Two of the Drosophila 
genes, however, have changed their function. These are 
fushi-tarazu (ftz) and Hox3. The latter gene evolved into 
bicoid (bed), zerknullt (zen), (also referred to as zl) and z2 
[28,29]. These genes have lost their homeotic function and 
now have new expression patterns. In Drosophila, ftz acts as 
a pair rule gene [30], and bed/zen-genes are involved in axis 
determination and formation of extraembryonic membranes 
[31]. However, in more basally branching arthropods, the 
expression patterns of Hox3 and ftz are consistent with ca- 
nonical Hox-like domains [32-36]. 

As Hox gene expression in arthropods is quite diverse 
(reviewed in e.g. [27]), it is difficult to reconstruct the ances- 
tral pattern of expression within the clade without reference 
to an outgroup, such as the onychophorans [37-39]. 

Here we report on the sequences and embryonic gene 
expression profiles of the complete set of the ten Hox 
genes in the onychophoran Euperipatoides kanangrensis. 
The new data contribute to our understanding of how 
the highly derived arthropod body plans have evolved 
from a rather uniform onychophoran-like ancestor. 
Based on the onychophoran data we reconstruct the an- 
cestral arthropod Hox gene profile as far as possible and 
use this in a comparative approach to detect phylogen- 
etic signal. We find support for Mandibulata and Tetra- 
conata, but not Myriochelata and Atelocerata. 

Material and methods 

Species husbandry, embryo treatment, in situ 
hybridization, nuclei staining and data documentation 

The collection, laboratory maintenance and dissection of 
embryos of Euperipatoides kanangrensis specimens are 
described in [40]. Whole mount in situ hybridization was 



performed as described in [10]. The developmental stage 
of all embryos was determined by analyzing embryos 
stained with the nuclear dye DAPI (4-6-Diamidino-2-phe- 
nylindole) and comparing to the stage table published by 
[41]. Embryos were analyzed under a Leica dissection 
microscope equipped with either an Axiocam (Zeiss) or a 
Leica DC100 digital camera. Brightness, contrast, and 
color values were adjusted in all images using the image 
processing software Adobe Photoshop CS2 (Version 9.0.1 
for Apple Macintosh). 

RT-PCR and gene cloning 

Total RNA was isolated from freshly-dissected embryos of 
E. kanangrensis via TRIZOL (Invitrogen). Poly- A RNA was 
isolated with the PolyATtract mRNA isolation system III 
(Promega) and used to produce cDNA using the Superscript 
first strand synthesis system for RT-PCR (Invitrogen). Short 
fragments of the Hox gene orthologs of E. kanangrensis Sex 
combs reduced (Ek-Scr), fushi-tarazu (Ekftz) and Abdom- 
inal-B (Ek-Abd-B) were isolated via RT-PCR with degener- 
ate primers. For that purpose mRNA was isolated and 
cDNA was synthesized from complete embryos represent- 
ing all stages from the 1-cell stage up to stage 21 [41]. A list 
of the primers used is available from the authors upon re- 
quest. Initial PCR fragments were elongated via 3 '-RACE 
(for Ek-Scr, Ekftz and Ek-Abd-B) or 5 '-RACE (for Ek-Antp) 
with the GENE RACER KIT (Invitrogen) to obtain suffi- 
ciently long fragments for subsequent in situ hybridization 
experiments. These sequences are available in GenBank 
under the accession numbers FR865437 (Ek-Scr), FR865438 
(Ek-ftz), FR865439 (Ek-Antp), and FR865440 (Ek-Abd-B). 
We also screened two independently prepared embryonic 
transcriptomes of E. kanangrensis for the presence of Hox 
genes and found a single copy of each of the expected Hox 
gene orthologs. The embryonic transcriptomes were made 
from comparable stages as used for the RT-PCR screen (1- 
cell stage to stage 21) [41]. We discovered the complete 
open reading frames of all onychophoran Hox genes. These 
sequences are available under accession numbers HE979835 
(Ek-lab), HE979836 (Ek-pb), HE979837 (Ek-Hox3), HE979838 
(Ek-Dfd), HE979839 (Ek-Scr), HE979840 (Ek-ftz), HE979841 
(Ek-Antp), HE979842 (Ek-Ubx), HE979843 (Ek-abd-A_splice 
variant I), HE979844 (Ek-abd-A_splice variant II), and 
HE979845 (Abd-B). 

All fragments were cloned into the PCRII vector 
(Invitrogen). Sequences of all isolated fragments were 
determined from both strands by means of Big Dye 
chemistry on an ABI3730XL analyser by a commercial 
sequencing service (Macrogen, Korea). 

Phylogenetic analysis 

The complete homeodomains were aligned in Clustal_X 
[42,43] and accuracy of the resulting alignment was 
checked by hand. Maximum likelihood analysis was 
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performed using the LG substitution model [44] as imple- 
mented in PhyLM [45]. 



Segmental nomenclature 

In order to facilitate description and comparison of the 
data, we label segments that express Hox genes in arthro- 
pods with the prefix "H"; HI is thus the most anterior seg- 
ment (corresponding variously to the onychophoran slime 
papilla segment, chelicerate pedipalps, insect and myria- 
pod "intercalary" and crustacean second antenna seg- 
ments) that expresses Hox genes (i.e. lab, pb and Hox3), 
thus avoiding the potentially thorny problem of how many 
segments lie in front of this. Onychophoran segments in 
general are labelled fap (frontal appendage), j (jaw), sp 



(slime papilla) and then LI- 15 for the walking limb- 
bearing segments. LI is thus H2. 

Results 

Sequence analysis of the ten onychophoran Hox genes 

Partial sequences of the E. kanangrensis Hox gene ortho- 
logs Ek-lab, Ek-pb, Ek-Hox3, Ek-Dfd, and Ek-Ubx have 
been described in [40,46]. We have now recovered the 
complete protein coding sequences of the ten onychoph- 
oran Hox genes. Orthology of the complete protein se- 
quences of the E. kanangrensis Hox genes was determined 
by comparison with published orthologs of the beetle Tri- 
bolium castaneum. The published sequence of ftz in the 
sea spider Endeis spinosa was added to the analysis to re- 
veal orthology of ftz genes (Additional file 1: Figure SI). 
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Figure 1 Expression of onychophoran Hox genes (lateral/ventral views). (A-H) Lateral views. (I-J) Ventral views. In all cases anterior is to the 
left. Upper row: schematic drawing of Hox gene expression. Middle row: bright-field photography (A-J). Lower row: DAPI staining (A'-J') of the 
same embryos as shown in exactly the same position as in the middle row (A-J). Darker shading in schematic drawings indicates stronger 
expression. Lighter shading indicates weaker expression. Note that the expression of Ek-Ubx mRNA forms a posterior to anterior gradient with 
unclear anterior expansion (marked with a question mark (?)). Arrows in l-J point to expression in the anal valves. Abbreviations: e, eye; j, jaw; 
LI -LI 5, first to fifteenth walking limb; saz, segment addition zone; sp, slime papilla. 
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We provide the alignment of conserved regions of the 
Hox genes from various arthropod and onychophoran 
species (see supplementary data for further information). 
Within the highly conserved regions lie diagnostic amino 
acids that are characteristic for each Hox gene (Additional 
file 2: Figure S2). We did not detect any duplicated Hox 
genes, neither in our PCR screens nor in the sequenced 
transcriptomes although two Abd-B paralogs have been 
reported for another onychophoran species, Akanthokara 
kaputensis (accession numbers AF011273 (Abd-B-1) and 
AF011274 (Abd-B-2)). Interestingly, we found one splice 
variant of Ek-abd-A with a hexapeptide (HX) sequence 
(£/c-<2&<i-A-variant I), and one without this highly- 
conserved motif (Ek-abd-A -variant II). For all other ony- 
chophoran Hox genes we only identified transcripts with 
the HX motif (except for Ek-Abd-B, which generally lacks 
a HX motif) (Additional file 2: Figure S2). This result dif- 
fers from the earlier published fragment of Ek-Dfd [40], 
which lacks a HX motif and also differs considerably in its 
complete N-terminal region from the newly recovered 
fragment of Ek-Dfd. The newly recovered fragment has 
significant sequence similarity with Dfd genes of other ar- 
thropods (Additional file 2: Figure S2). The previously 
published short fragment of Akanthokara Dfd [47] most 
probably represents an insect sequence since it contains a 
number of insect- specific amino acids (Additional file 2: 
Figure S2). The earlier published sequence of Ek-Hox3 



[40] also differs from the newly-recovered fragment. The 
former sequence contains a short string of additional 
amino acids between the HX motif and the homeodomain 
(HD). This sequence was found neither in the sequenced 
transcriptomes nor in newly-cloned fragments recovered 
by means of RT-PCR with gene specific primers that 
flanked the sequence in question (not shown). It remains 
unclear, however, whether the former sequence represents 
an artefact or a rare splice variant. 

Expression of the ten onychophoran Hox genes 

The expression profiles of Ek-lab, Ek-pb, Ek-Hox3 and Ek- 
Dfd have been described in [40] (Figures 1A-D and 2A-D). 
The expressions of all ten Hox genes are continuous from 
their respective anterior borders of expression back to the 
rear of the developing embryos. 

labial 

In addition to the previously-described expression we find 
that Ek-lab is expressed considerably more strongly in the 
slime papillae-bearing (sp) segment (Figures 1A and 2A) 
as compared to other segments, and expression is up- 
regulated in certain regions of the developing neuroecto- 
derm of the trunk (Figure 2A). A spot of Ek-lab expression 
is located anteriorly and proximally in the slime papillae 
(Figure 2A). 
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Figure 2 Expression of onychophoran Hox genes (ventral views). In all cases anterior is to the left. Arrows in A and B point to up-regulated 
expression of lab and pb in the ventral nervous system. The asterisk (*) in A indicates damage of ventral tissue as a result of the hybridization 
procedure. See text for further information. Abbreviations; av, anal valves; fap, frontal appendage; hi, head lobe; j, jaw, LI -LI 5, first to fifteenth 
walking limb; sp, lime papilla, vee, ventral extraembryonic ectoderm. 
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Figure 3 Anterior borders of Ek-Hox3, Ek-Scr, Ek-Dfd and Ek-ftz 
expression. A) Ek-Hox3 expression in sp and LI. Asterisk marks 
anterior border of expression in the middle of the slime papilla (sp) 
segment (cf. Eriksson et al. 2010). B) Expression of Ek-Dfd extends to 
the anterior of the L1 segment (cf. Eriksson et al. 2010). C) Ek-Scr. 
Asterisk marks anterior border of expression in the middle of the L1 
segment. D) Ek-ftz. Asterisk marks anterior border of expression in 
the middle of the L2 segment. Abbreviations as in Figure 1. Dashed 
lines indicate segmental borders. 



proboscipedia 

The only difference between the published expression 
pattern and our new data is that Ek-pb is up-regulated 
in some regions of the neuroectoderm (Figure 2B). 

Hox3 

Ek-Hox3 remains expressed in the posterior half of the 
slime papillae, and does not later disappear from this tis- 
sue (Figure 2C) cf. [40]. 

Deformed 

The anterior border of Ek-Dfd is at the border between 
the jaw-bearing segment and the first trunk segment (LI) 
(Figures ID and 2D). In embryos younger than stage 11, ex- 
pression only extends into L2 (Additional file 3: Figure S3). 

Sex combs reduced 

The anterior border of Ek-Scr expression lies approxi- 
mately in the middle of the first walking limb-bearing 
segment (LI) (Figures IE, 2E and 3). From there the ex- 
pression extends throughout the complete posterior 
body including the anal valves. 

fushi-tarazu 

The anterior border of Ek-ftz lies in the middle of the sec- 
ond trunk segment (L2) (Figures IF, 2F and 3). In embryos 
younger than stage 12, expression only extends anteriorly 
into L3 (Additional file 3: Figure S3). 

Antennapedia 

Ek-Antp is expressed in all posterior segments including all 
of L4. Expression in L4, however, is considerably weaker than 
in more posterior segments. In segments L5-L7, Ek-Antp is 
expressed stronger than in L4, but still slightly weaker than 
in L8 to L15, which is the last segment (Figures 1G and 2G). 
Whether expression in L5 to L7 forms a gradient, or if it is at 
the same level in these segments, is not revealed by our in 
situ hybridization technique. As described for Scr and ftz, ex- 
pression of Antp also shifts towards the anterior at around 
stage 13 (Additional file 3: Figure S3). 

Ultrabithorax 

The Ek-Ubx gene is clearly expressed in an anterior to 
posterior gradient. The most anterior segment expressing 
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Ek-Ubx at a detectable level is L6 (Figures 1H and 2H). 
The gradient makes it difficult to determine unambigu- 
ously the anterior- most extent of Ek-Ubx expression. 
Low-level transcription (below detectable range) may be 
present in L5. 

abdominal-A 

Ek-abd-A is exclusively expressed in the last trunk seg- 
ment (L15) and the anal valves (Figure II). Notably, ex- 
pression in the mesoderm is stronger than in the 
ectoderm and at late developmental stages Ek-abd-A is 
strongly expressed in the interface between L15 and the 
anal valves (Additional file 4: Figure S4). Note the differ- 
ence between this specific staining and frequently occur- 
ring unspecific signal in L15 (Additional file 5: Figure S5). 

Abdominal-B 

Ek-Abd-B is expressed in the mesoderm of the entire anal 
valves, but ectodermal expression is restricted to the very 
posterior of the anal valves (Figure 1J and Additional file 
4: Figure S4). 

Generally, none of the ten onychophoran Hox genes are 
expressed in tissue dorsal to the base of the appendages 
(Figure 1). This is unlike the situation in arthropods where 



expression of the Hox genes is generally not restricted to 
ventral tissue e.g. [27,34,48]. 



Discussion 

Comparison of arthropod and onychophoran Hox gene 
expression data 

Although the expression patterns of the Hox genes are 
relatively conserved across the arthropods, some clade- 
specific differences exist [27]. Comparison between arthro- 
pod and onychophoran Hox gene expression patterns may 
aid in reconstructing ancestral and derived arthropod ex- 
pression: firstly, shared features of onychophorans and ar- 
thropods are likely to have been present in the last 
common ancestor of both; and secondly, the onychoph- 
oran data can in some instances be used to polarise the 
order of Hox gene change within the arthropod tree itself. 
We have therefore mapped Hox gene expression patterns 
onto an arthropod plus onychophoran tree (Figure 4). The 
Mandibulata grouping (myriapods plus insects plus crusta- 
ceans) seems favoured over Paradoxopoda (myriapods and 
chelicerates) in recent phylogenies [49-53] and thus we 
have chosen this topology to map our data onto, although 
as noted below, some aspects of the onychophoran data 
do support Mandibulata and none support Paradoxopoda. 
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Figure 4 Summary of phylogenetic signals based on onychophoran and arthropod Hox gene expression. Conserved and diverged characters 
are plotted on the nowadays best-supported arthropod phylogeny. Onychophora, and onychophoran data, serve as outgroup. Dark red boxes denote 
events during arthropod Hox gene evolution. Light red boxes denote unresolved character-state; in these cases the character-state cannot be 
determined ancestral or derived because outgroup data for a clade Onychophora + Arthropoda are not available. 
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Figure 5 Comparison of Hox gene expression profiles in the 
onychophoran E. kanangrensis, the myriapod Glomeris marginata 
and the harvestman Phalangium opilio, and reconstruction of the 
putative Hox gene pattern in the common ancestor. Orthologous 
genes are with identical colours. Conserved anterior and posterior 
borders of expression (compared to E konongrensis) are marked with 
asterisks (*). Borders of expression that are conserved in the myriapod 
and the harvestman (but not the onychophoran) are marked with 
filled circles (•). Smaller bars indicate weaker expression. 'Delayed' 
expression is marked with 'd]'. Green bar to the left indicates good 
support for the predicted corresponding expression patterns in the 
ancestor; yellow bars indicate ambiguous support for the 
reconstructed expression patterns; red bar indicates that the ancestral 
expression pattern of abd-A is not predictable. This is also indicated by 
question marks (?). Dashed vertical lines indicate an either 
segmental or parasegmental (dashed line) anterior expression 
border. Hox gene-expressing segments are labelled with the prefix 
'H'. Abbreviations are as in Figure 1, and an, antennal segment; ch, 
cheliceral segment; L1-L4 (in chelicerates), walking limb-bearing 
segments one to four; md, mandibular segment; mx, maxillary 
segment; 01-09, first to ninth opisthosomal segment; pa, pre-antennal 
region; ped, pedipalpal segment; pmd, premandibular segment; pmx, 
postmaxillary segment; SAZ, segment addition zone; T1-T8, first to 
eighth walking limb-bearing segment. 



Onychophorans and arthropods share several import- 
ant aspects of their Hox gene expression patterns. 
Firstly, in both onychophorans and arthropods, Hox 
gene expression is absent from tissue anterior to the 
slime papillae-bearing (intercalary/premandibulary/ped- 
ipalpal) segment [25,34,40,54]: i.e., Hox gene expression 
never extends to the anterior-most segments in either 
clade. Secondly, all Hox genes, except onychophoran 
Ubx, have a distinct anterior expression border, al- 
though at early developmental stages in members of 
both phyla this border may be located one segment 
more posteriorly. 

A difference between onychophorans and arthropods is 
that in arthropods Scr is expressed in H3, but in ony- 
chophorans it is in H2. It would appear that either the 
onychophoran pattern is derived, or that a shift in the an- 
terior extension of Scr evolved in the lineage leading to 
the arthropods (Figures 4 and 5 and Additional file 6: 
Table SI). Another apparent difference concerns the pos- 
terior extent of gene expression. In E. kanangrensis the ex- 
pression of all Hox genes extends back to the rear of the 
embryo [40], this study; but in arthropods, the "anterior" 
Hox genes are typically not expressed in the posterior part 
of the body. However, data on Hox gene expression in 
species of the Chelicerata and Myriapoda provide evi- 
dence that the onychophoran Hox gene pattern may in 
this respect reflect the ancestral condition for Onychoph- 
ora + Arthropoda. In the harvestman Phalangium opilio 
[36] (Figure 5A), the spider Achaearanea tepidariorum 
[55], the mite Archegozetes longisetosus [56], the centipede 
Lithobius atkinsoni [54] and the millipede Glomeris mar- 
ginata [34] (Figure 5A), the expression of at least some 
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anterior Hox genes extends from their anterior border 
of expression to the rear end of the embryo. In other 
cases anterior Hox genes are also expressed in the most 
posterior segment(s) and/or the segment addition zone 
(SAZ) (Figure 5A). We believe that weak (but clearly de- 
tectable) expression of some Hox genes in the posterior 
segments and expression in the last segment (or the 
SAZ) may represent evolutionary remnants of an ances- 
tral state more clearly retained in onychophorans. Taken 
together, these data indicate a tendency in arthropod 
evolution to restrict anterior Hox gene expression to 
distinct anterior regions. 

Clade-specific comparisons between onychophorans 
and arthropods 

As well as the general similarities (and differences) 
highlighted above, onychophorans share several similarities at 
the arthropod clade-specific level, in particular in conserva- 
tion of anterior Hox gene expression borders (summarised in 
Figure 4). Although the anterior Hox gene pattern of pb and 
Hox3 is highly conserved, it deviates in some clades at least 
in insects, and these latter seem therefore to be derived. Sec- 
ondly, deviation from the pattern of ftz expression found in 
onychophorans is only found in crustaceans and insects, sup- 
porting the Tetraconata (Pancrustacea) hypothesis. Thirdly, 
the anterior border of Antp seen in onychophorans has 
switched from H5 to H4 in myriapods, crustaceans and 
insects, in agreement with the Mandibulata concept. 
Fourthly, the posterior reduction of expression of the 
anterior Hox genes discussed above progresses from ony- 
chophorans, where expression of all the Hox genes extend 
to the rear end, to crustaceans and insects, where ex- 
pression is more confined to distinct anterior regions 
(Figure 5 A and Figure twelve in [27]). This again is in 
agreement with the Tetraconata concept. Conversely, 
the polarisation possible with the onychophoran data 
does not support either Myriochelata or Atelocerata. 

Taken together, the differences and similarities in the ex- 
pression data between onychophorans and arthropods allow 
a reasonable degree of reconstruction of the ancestral Hox 
expression for the entire clade (Figure 5B). We find good 
support for the reconstruction of the pattern of anterior Hox 
genes {lab, pb, Hox3, Dfd, Scr, Antp), but also that the recon- 
struction of posterior Hox genes is difficult (Ubx, abd-A, 
Abd-B) (Figure 5B). The varying anterior border of abd-A in 
the onychophoran and various arthropods makes it impos- 
sible to reconstruct its ancestral pattern. The lability of abd-A 
is highlighted by the loss of this gene from e.g. the mite Tet- 
ranychus urticae [57]. Among the reconstructable patterns 
some uncertainties remain, such as whether the Ubx gradient 
in onychophorans is ancestral or derived, and the anterior ex- 
pression boundary of Scr. These uncertainties might in 
principle be resolved by reference to a suitable outgroup such 
as a tardigrade or a cycloneuralian worm. Unfortunately, Hox 



gene expression and action in the nematode Caenorhabditis 
elegans is highly derived [58], and no other cycloneuralian ex- 
pression patterns are known. Development of a priapulid in 
situ hybridization protocol [59], however, opens the possibil- 
ity of such data being obtained in the future. 

Hox gene expression and tagmosis in arthropods 

Hox gene expression has long been associated with the 
evolution and maintenance of tagmosis, i.e. the character- 
istic grouping or fusion of functionally similar segments 
e.g. [24,25,36,60]. Ever since the earliest days of Hox gene 
research, the dominant hypothesis has been that the Hox 
genes primarily act to specify regions of the body, includ- 
ing above all the tagmata [16]. 

Although onychophorans have a relatively unspecialised 
body plan, with only the head and most posterior segment 
being differentiated, their Hox gene expression patterns 
show a surprising degree of sophistication, for example in 
the gradient of expression of Ubx. This pattern mirrors a 
similar one documented in limb specification, which in on- 
ychophorans is also surprisingly sophisticated [10]. 

In the onychophoran studied, the anterior segments (sp- 
L5) and the last segment (LI 5) express a unique set of 
Hox genes, which allows for the characterization of each 
of these segments sensu [17] (Figure 5 A). The posterior 
segments (L6-L14), however, all express the same eight 
Hox genes. The same is true for many arthropods where 
the posterior segments are usually not characterised by a 
unique set of Hox genes [27]. It is possible that posterior 
segments in onychophorans do not require a specific set 
of Hox gene input as they appear to be morphologically 
identical, but that would then raise the question why the 
anterior (morphologically identical) trunk segments in 
onychophorans do express distinct sets of Hox genes 
(Figure 5A). It is of course possible that there are some 
cryptic morphological differences in these segments that 
are being regulated by differential Hox gene expression. 
However, if this is the case, then it is possible that such 
cryptic morphological differences are present in the pos- 
terior segments too, and this would raise the question of 
how these segments are specified within identical Hox 
gene expression domains. Generally, the mRNA expres- 
sion pattern can differ from the protein pattern as a result 
of translational repression, and this could lead to different 
Hox-protein and Hox-mRNA landscapes. Further differ- 
ences in segmental Hox gene patterning might also in- 
volve Hox cofactors [61-63], the presence of alternative 
splice variants e.g. [64] (possibly as represented by the 
different splice variants of Ek- abd-A), and temporal dif- 
ferences in Hox gene expression [65,66]. Such temporal 
differences include the shifting anterior expression bor- 
ders of Ek-Scr, Ek-ftz, and Ek-Antp (Additional file 3: 
Figure S3). 
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Ek-Ubx mRNA is expressed in an anterior to posterior 
gradient with detectable transcripts from the very pos- 
terior of the embryo to at least L6. It is thus possible 
that the different levels of Ek-Ubx mRNA (and resulting 
£/c-Ubx protein) are sufficient to give these segments a 
unique Hox signature. Studies in at least Drosophila do 
indeed show that Ubx function in segmentation and 
limb development is dependent on different levels of ex- 
pression [67,68]. Similarly, in the crustacean Parky ale, 
different levels of Ubx expression may be responsible for 
the development of different types of appendages, and 
knock down of Ubx results in homeotic transformation 
of these segments [69]. The Ek-Ubx gradient thus offers 
one possibility of how each posterior segment is 
uniquely patterned, even though they express the same 
set of Hox genes. 

Additional files 



the ectoderm of the anal valves. Abbreviations: a, anus; av, anal valve; ec, 
ectoderm; mes, mesoderm; L14 and L15, fourteenth and fifteenth 
walking-limb bearing segment. 

Additional file 5: Figure S5. Frequently occurring staining artefacts. In 
all panels anterior is to the left. A) Stage 19 embryo. Ventral view. Typical 
staining artefacts are often observed in older stage embryos where the 
surface of the frontal appendages attracts unspecific signal. A similar false 
positive signal appears regularly in the opening of the slime papillae. The 
coelomic cavity of the L15-segment often stains unspecifically. Long 
staining time frequently results in unspecific staining of the yolk. B/C) A 
typical artefact: staining of one of the two coelomic cavities (arrowheads) 
of the head in a stage 13 and a stage 20 embryo. Abbreviations: fap, 
frontal appendage; L15, fifteenth walking limb; sp, lime papilla. 

Additional file 6: Table SI. Summary of available literature on Hox 
gene expression in arthropods and onychophorans. Compared is the 
anterior-most extension of Hox gene mRNA (in some cases protein) 
expression. 
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